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' Abstract 

, Some time ago, Thies et al. showed that the Gross-Neveu model with a bare mass term 

possesses a kink-antikink crystalhne phase. Corresponding self-consistent solutions, known 
earlier in polymer physics, are described by a self-isospectral pair of one-gap periodic Lame 
potentials with a Darboux displacement depending on the bare mass. We study an unusual 
supersymmetry of such a second order Lame system, and show that the associated first order 
Bogoliubov-de Gennes Hamiltonian possesses the own nonlinear supersymmetry. The Witten 
index is ascertained to be zero for both of the related exotic supersymmetric structures, each 
' of which admits several alternatives for the choice of a grading operator. A restoration of the 

<~i , discrete chiral symmetry at zero value of the bare mass, when the kink-antikink crystalline 

condensate transforms into the kink crystal, is shown to be accompanied by structural changes 
, in both of the supersymmetries. We find that the infinite period limit may or may not change 

^ ' the index. We also explain the origin of the Darboux dressing phenomenon recently observed 

0^ , in a non-periodic self-isospectral one-gap Poschl- Teller system, which describes the Dashen, 

(N . 

Hasslacher and Neveu kink-antikink baryons. 

in 

(N 

1 Introduction 



O 



Oh: 



The Gross-Neveu (GN) model [UEIIH] is a remarkable (l-l-l)-dimensional theory of self-interacting 
^ ' fermions that has no gauge fields or gauge symmetries, but exhibits some important features of 

quantum chromodynamics, namely, asymptotic freedom, dynamical mass generation, and chiral 
symmetry breaking [3]. It has been widely studied over the years and the richness of its properties 
is still astonishing. Some time ago, Thies et al. showed that at finite density, the ground state of 
the model with a discrete chiral symmetry is a kink crystal [5], while the kink-antikink crystalline 
phase was found in the GN model with a bare mass term p]. Then, Dunne and Basar derived 
a new self-consistent inhomogeneous condensate, the twisted kink crystal in the GN model with 
continuous chiral symmetry [Tj |8]. On the other hand, the relation of the GN model with the 
sinh-Gordon equation and classical string solutions in AdSa has been observed recently [9l llOj. 

These two classes of the results seem to be different, but both are rooted in the integrability 
features of the GN model, and may be related to the Bogoliubov-de Gennes (BdG) equations 
incorporated implicitly in its structure. It is because of these properties that the model finds 
many applications in diverse areas of physics. Particularly, the model has provided very fruitful 
links between particle and condensed matter physics, see [HI [12] and [13]. 
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The origin of the model itself may also be somewhat related to the BdG equations. We briefly 
discuss these equations to formulate the aim of the present paper. 

The BdG equations [14j in the Andreev approximation [15] is a set of two coupled linear dif- 
ferential equations, which can be presented in a form of a stationary Dirac-type matrix equation, 

Gill) = ojip, Gi=aai-- a2^{x) . (l-l) 

I ax 

The scalar field A(a;) is determined via a self-consistency condition, which often referred to as 
a gap equation. Equation (jl.ip arose in the theory of superconductivity by linearizing the non- 
relativistic energy dispersion (in absence of magnetic field), or, equivalently, by neglecting second 
derivatives of the Bogoliubov amplitudes, see |16] . A constant a is proportional there to the 
Fermi momentum hkp. In what follows we put a = 1 and h= 1. 

The Lagrangian of the GN model of the A'^ species of self-interacting fermions is 

CGN = H^l''^^-mo)i; + ^g\^P^P)\ (1.2) 

where g"^ is a coupling constant, the summation in the flavor index is suppressed, and a bare 
mass term ~ mo, which breaks explicitly the discrete chiral symmetry ip — t- j^ip of the massless 
model, is included 0. It is the two-dimensional version of the Nambu-Jona-Lasinio model [T7] 
[with continuous chiral symmetry reduced to the discrete one] . The latter is based on an analogy 
with superconductivity, and was introduced as a model of symmetry breaking in particle physics. 
There are two equivalent methods to seek for solutions for the GN model. One of them is the 
Hartree-Fock approach, in which self-consistent solutions to the Dirac equation {i'y^d^ — S)^l^ = 
are looked for, with spinor and scalar fields subject to a constraint of the form {S{x) — m-o) = 
—Ng'^{xpip), see [UIUITI]. For static solutions, under appropriate choice of the gamma matrices, 
the Dirac equation takes a form of the BdG matrix equation (jl.ip . with Gi as a single particle 
fermionic Hamiltonian. The condensate field S{x) is identified with a gap function A(x), while 
the constraint corresponds to the above mentioned gap equation. Another approach to seek 
solutions for the GN model, in which the BdG equations also play a key role, is via a functional 
gap equation |19ll20j . There, the condensate field is given by stationary points of effective action, 
and a connection of the GN model with integrable hierarchies can be revealed, see [3 El [20l 121] . 
In light of this, the relation of the GN model to the sinh-Gordon equation does not seem to be 
so surprising as the BdG equations arise (in a slightly modified form) as an important ingredient 
in solving the sine-Gordon equation, see [22ll23] . 

We now return to the BdG matrix system (jl.ip . By squaring, the equations decouple, 

(P 

Hi> = Ei>, E = u;'^, H = --^ + A'^ - asA' . (1.3) 

From the viewpoint of the second order system H = G\, the first order matrix operator Gi is a 
nontrivial integral of motion, [H,Gi\ = 0. Having also an integral (T3, [i^,(73] = 0, which anti- 
commutes with Gi , we obtain a pattern of supersymmetric quantum mechanics with (T3 identified 
as a grading operator. Though a system of the first and second order equations (jl.ip and (jl.3p was 
exploited in investigations on superconductivity, its superalgebraic structure, which also includes 
the second supercharge G2 = ia^Gi, seems to have gone unnoticed before the theoretical discovery 

^The investigation of model (|1.2p is motivated in [6] by a massive nature of quarks; there, the 't Hooft hmit 

A'^ — Qo, Ng^ — const, is considered. 
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of supersymmetry in particle physics. Supersymmetric quantum mechanics was then developed 
by Witten as a toy model for studying the supersymmetry breaking in quantum field theories 
|24j . Later, the relation of supersymmetric quantum mechanics with Darboux transformations 
was noticed [25], and found many applications [26]. 

Braden and Macfarlane [27], and, in a broader context, Dunne and Feinberg [28] observed 
that the Darboux transformed, supersymmetric partner of the one-gap periodic Lame system 
|29j with a zero energy ground state is described by the same potential but translated for a half- 
period. The superpartner, therefore, also has a zero ground state. Such a system is described by 
unbroken supersymmetry, in which, however, the Witten index takes zero value. For a class of su- 
perpesymmetric systems with super-partner potentials of the same form a term self-isospectrality 
was coined by Dunne and Feinberg [28] . The supersymmetric Lame system considered in [27^ [2H] 
corresponds to the kink crystalline phase discussed in [5j , which describes a periodic generalization 
of the Callan- Coleman- Gross- Zee (CCGZ) kink configurantions of the GN model, see [S] [TSl [30] 
and [16] , It was known earlier as a self-consistent solution to the GN model in the context of 
condensed matter physics [3Tj , see also [32l [33l IS] . 

The Lame system, like non-periodic reflectionless solutions of the GN model, belongs to a 
special class of the finite-gap systems [25l [35] El- Some time ago, it was found that such systems 
in an unextended case (i. e. when a second order Hamiltonian has a single component), are 
characterized by a hidden, peculiar nonlinear supersymmetry [ 37[ 138] . It is associated with a 
corresponding Lax operator (integral), and the grading is provided there by a reflection operator. 
As a consequence, supersymmetric structure of an extended system [with a matrix Hamiltonian 
of the form (II. 3p ] turns out to be much richer than that associated with only the first order 
supercharges Ga, a= 1,2, and integral a^, see [39]. It has also been shown recently [40j that the 
self-isospectral Poschl- Teller system (PT), which describes the Dashen-Hasslacher-Neveu (DHN) 
kink-antikink baryons [2J, is characterized by a very unusual nonlinear supersymmetric structure 
that admits six more alternatives for the grading operator in addition to the usual choice of 
(73 . All the local and non-local supersymmetry generators turn out to be the Darboux-dressed 
integrals of a free non-relativistic particle. Moreover, it was shown there that the associated BdG 
system, with the matrix operator (jl.ip identified as a first order (Dirac) Hamiltonian, possesses 
its own, nontrivial nonlinear supersymmetry. 

In the present paper we investigate the exotic supersymmetric structure of the kink-antikink 
crystal of [6^ i31j; which is a self-consistent solution of the GN model (|1.2p with a real gap 
function A(x; r). Parameter r is related to niQ and controls a central gap in the spectrum of the 
first order BdG Hamiltonian operator (jl.ip . Simultaneously, it defines a mutual displacement, 
2r, of superpartner Lame potentials in correspondence with the structure of the second order 
Schrodinger operator (jl.Sp . One more parameter, not shown explicitly here, defines a period of 
the crystal. A quarter-period value of r corresponds to the kink crystal solution of [5] for the 
model ()1.2p with niQ = 0, which was considered in [27[ 128) . We also study different forms of the 
infinite-period limit applied to the supersymmetric structure. A priori the picture of such a limit 
has to be rather involved : the Darboux dressing relates the non-periodic kink-antikink system 
to a free particle, while the Darboux transformations in the periodic case are expected to be just 
self-isospectral displacements, see [3T ] [39 ] W\ ] [^. 

The outline of the paper is as follows. In the next section, we discuss the main properties of the 
one-gap Lame system. In section 3 we construct its self-isospectral extension by employing certain 
eigenfunctions of the Lame Hamiltonian. We investigate the action of the first order Darboux 

^There is also a relation of one-gap Lame equation with the sine-Gordon equation, see |36) . 



3 



displacement generators, and discuss the spectral peculiarities of the obtained supersymmetric 
system. Section 4 is devoted to the study of the properties of a superpotential (gap function) that 
is an elliptic function both in a variable and a shift parameter. These properties are employed 
in section 5, where we construct the second order intertwining operators, identify further local 
matrix integrals of motion, and compute a corresponding nonlinear superalgebra. In section 6 
we show that the system possesses six more, nonlocal integrals of motion, each of which may be 
chosen as a Z2 grading operator instead of a usual integral (T3 of the supersymmetric quantum 
mechanics. We discuss alternative forms of the superalgebra associated with these additional 
integrals and their action on the physical states of the system. In section 7, we investigate a 
peculiar nonlinear supersymmetry of the associated first order BdG system. Section 8 is devoted 
to the infinite period limit of the both, second and first order supersymmetric systems. In section 
9 we clarify the origin of the Darboux dressing phenomenon that takes place in the non-periodic 
self-isospectral PT system, that was revealed in [30]. In section 10 we discuss the obtained results. 
To provide a self-contained presentation, the necessary properties of Jacobi elliptic functions and 
of some related non-elliptic functions are summarized in the two appendices. 



2 One-gap Lame equation 

In this section we discuss the properties of the Lame system which is necessary for further 
constructions and analysis. 

Consider the simplest (and unique) one-gap periodic second order system described by the 
Lame Hamiltonian 

H = + 2k^sn'^x - k'^ . (2.1) 

An additive constant term is chosen here such that a minimal energy value (the lower edge of 
the valence band, see below) is zero. Potential V{x) = 2A;^sn^x — A:-^ is a periodic function with 
a real period 2K (and a pure imaginary period 2iK') H. The general solution of the equation 

H^{x) = E^{x) (2.2) 

is given by [29] 

*±(^) = ^%^exp[TxZ(a)] . (2.3) 

Here H, and Z are Jacobi's Eta, Theta and Zeta functions, and the eigenvalue E = E[a) is 
defined by the relation 

E{a)=di?a. (2.4) 



The Hamiltonian (j2.ip is Hermitian, and we treat (j2.2p as the stationary Schrodinger equation 
on a real line. We are interested in the values of the parameter q, which give real E. dn^a is 
an elliptic function with periods 2K and 2iK', and its period parallelogram in a complex plane 
is a rectangle with vertices in 0, 2K, 2K + 2iK' and 2zK'. We then look for those a in the 
period parallelogram for which dna takes real or pure imaginary values. They can be taken, for 
instance, on the border of the rectangle shown on Fig. [TJ We have, particularly, 

E{V. + il3) = k'^cu^{l3\k')iid^{l3\k'), 0</3<K', A;'^ > S(K + i/3) > , (2.5) 

E{il3)=du^{P\k')iic^{p\k') = k'^ + k'^iic^{l3\k'), 0</3<K', I < E{ip) < 00 . (2.6) 



See Appendices A and B for notations and properties we use for Jacobi elliptic and related functions. 
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zK' J -oo<E<0 K+iK' 



0<E<k'^ 



0' k'^<E<l 

Figure 1: The sides of the rectangle are mapped by (j2.4p onto the indicated energy intervals. 
The vertical (horizontal) sides shown in green (red) correspond to the two allowed (forbidden) 
bands. Vertices a = K + iK', K and are mapped, respectively, into the edges E = 0, k''^, and 1 
of the valence, < E < /c'^, and conduction, 1 < E < oo, bands, which are described by periodic, 
dnx {E = 0), and antiperiodic, cnx (E = k'"^) and snx (i? = 1), functions. Vertex iK! as a limit 
point on a horizontal (vertical) side corresponds to E = —oo [E = +oo). 



1<E<^ 



For (|2.5|) and (|2.6|) . eigenfunctions in (j2.2p are bounded on a real line, that corresponds to the two 
allowed (valence and conduction) bands in the spectrum. In contrast, for a = (3 and a = /3 + zK', 
j3 G (0,K), a real part of Z(a) is nonzero, and eigenfunctions (12. 3p are not bounded for — )■ oo. 
This corresponds to the two forbidden zones, — oo < E <Q and k''^ < E < 1. 
Differentiation of (j2.5p and (j2.6p in (3 gives a relation 



— = 2v{E) , P{E) = E{E - k"){E - 1) . (2.7) 

The third order polynomial P{E) takes positive values inside the allowed bands, and turns into 
zero at their edges. r]{E) takes values —1 and +1 in the valence and conduction bands, respec- 
tively. 

Inside the two allowed bands, ()2.3p are quasi-periodic Bloch wave functions, 

^l{x + 2K) = e^*2K«{i?)^a ^ = ^ " ' (2-^) 

where a first term in quasimomentum (crystal momentum) n{E) originates from the imparity of 
H function. In the valence, (|2.5p . and conduction, (|2.6p . bands its values are given by 

k{E{V. + i/3)) = [Z{/3\k') + 21^/3 - A:'2cn(/3|A;')sn(/3|A:')nd(/3|A:')] , (2-9) 

K(i?(2/3)) = [Z{/3\k') + 21^/3 - dn(/3|fe')sn(/3|fe')nc(/3|A;')] . (2.10) 



With the help of (j2.4p and (j2.7p one finds a differential dispersion relation 

dK ,^,E-(E/K) 

— = 7]{E) , , 2.11 

dE ^ 2/POB) 

where E is a complete elliptic integral of the second kind, see (IB.ip . Taking into account a relation 
k'^ < ^ < 1, see Appendix B, one finds that within the both allowed bands quasimomentum is 
increasing function of energy. It takes values and 7r/2K at the edges E = and E = k'"^ of the 
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valence band, where Bloch-Floquet functions reduce to the periodic, dnx, and antiperiodic, cnx, 
functions in the real period 2K of the system. Within the conduction band, quasi-momentum 
increases from 7r/2K to +00. At the lower edge E = 1, two functions ()2.3p reduce to the 
antiperiodic function snx. At all three edges of the allowed bands, derivative of quasimomentum 
in energy is +00. For large values of energy, E — )• +00, we find that k{E) ~ ^/E, i.e. Bloch 
functions ()2.3p behave as the plane waves, 

Second, linear independent solutions at the edges of the allowed bands Ei = 0, k''^, 1 are 
^j(x) = 'ipi{x)Ii, where !{ = J dx / ipf {x) , and tpi = dnx, cnx, snx, i = 1,2,3. The integrals 
are expressed in terms of non-periodic incomplete elliptic integral of the second kind ()B.2p . 
Xi = ;^E(x + K), X2 = 2; - ^E(x + K + iK'), X3 = x - E(x + iK'). ^i(x) are not bounded on 
M and correspond to non-physical states. These non-physical solutions follow also from general 
solutions (|2.3p . For instance, ^3(x) may be obtained as a limit of (^'"(x) — ^'°(x))/a as a — )• 0. 
Eq. (12. 3p provides a complete set of solutions for (12. 2p as the second order differential equation. 
Notice also that Bloch states ()2.3p within the allowed bands are related under complex conjugation 
as (^'"(x))* = T/^'"(x), where r] is the same as in ()2.7p . 

In conclusion of this section we note that the function P{E) in Eqs. (|2.7p . (|2.1ip is a spectral 
polynomial. It will play a fundamental role in a nonlinear supersymmetry we will discuss below. 



3 Self-isospectral Lame system 

Consider the lower in energy E forbidden band by extending it with the edge value = of 
the valence band. We introduce a notation — 2r -|- iK' for the parameter a that corresponds to 
the extended interval —00 < E <{). By taking into account relations dn {—u) = dn (u -|- 2K) = 
— dn (ii-|-2iK') = dnu, it will be convenient do not restrict the values of r to the interval [— K/2, 0), 
but assume that r € M, while keeping in mind that E — )• — c« for r — t- nK, n G Z. After a shift 
of the argument x — )• x-|-t, the corresponding function \I'" from (j2.3p with a = — 2r-|-iK' takes, 
up to an inessential multiplicative constant, the form 

||-jexp[xz(r)] = F(x;r), (3.1) 

where we have introduced the notations x+ = x -|- r, x_ = x — r, 

z(r) = -iK(E(-2T + iK')) = ?(r) + Z(2r) = --^ ln(e(2T) sn2r) , (3.2) 

2 dr 

qir) = — — — lnsn2r = ns2r cn2Tdn2r . (3.3) 

^ ^ 2dT ^ ^ 

F{x; r) is a quasi-periodic in x and periodic in the r function, F{x + 2K; r) = 
exp (2Kz(r))F(x; r), F{x; r-|-2K) = F{x; r). It is a regular function of r, save for r = nK, n G Z, 
[which correspond to the poles a = 2nK -|- iK' of dna in (|2.4p ]. where F{x; r) with x 7^ under- 
goes infinite jumps from to +00. Since z(K/2) = 0, function (13. ip reduces at r = K/2 (up to an 
inessential multiplicative constant) to a periodic in the x function dn (x -|- ^K) which describes a 
physical state with energy E = at the lower edge of the valence band of the system H{x + ^K). 
F{x;t) is a nodeless function that obeys the relations F{x; — r) = F{—x;t) = 1/F{x;t) and 

[H{x+)+e{T)]F{x;T) = 0, where e(r) = -^(-2r iK') = cn^ 2t ns^ 2r . (3.4) 
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Define a first order differential operator 

V{x; r) = F{x; r) ^ = ^ - A(x; r) , V\x; r) = -V{x; -r) , (3.5) 

where 

A(x;r) = |;^. (3.6) 

Operator ()3.5p annifiilates function (j3.ip . 2?(a;; r)-F(x; r) = 0, and we find that 

V\x;T)V{x-T)=H{x+) + e{T), V{x-T)V\x]T) = H{x^) + e{T). (3.7) 

By virtue ofe(iK) = 0, a non-shifted Lame Hamiltonian operator ()2.ip factorizes then as -ff(x) = 
V [x + ^K; ^K) (x + ^K; ^K). The alternative product produces a shifted in the half-period 
K system, H{x + K) = [x + ^K) V {x + ^K; ^K) . It is this factorization of a pair of 
Lame Hamiltonians H{x) and -ff(x-l-K) that underlies a usual supersymmetric structure studied 
in [28] in the light of a phenomenon of self-isospectrality. 

Notice that while F(x_;t) is, up to a multiplicative constant, a non-physical eigenfunction 
^-2t+iK q£ i/(x) of energy — e(r), function -F(x+; — r) = l/F(x+;r) coincides, up to a mul- 
tiplicative constant, with another eigenfunction ^~2r+iX ^^-^ q£ ^^^j) with the same eigenvalue. 

According to (j3.7p . the mutually shifted Hamiltonians H{x + r) and H{x — r) form a super- 
symmetric, self-isospectral periodic one-gap Lame system 

n = diag {H{x+),H{x^)), (3.8) 

see Fig. [2l for which A(x; r) plays the role of the superpotential that obeys the Ricatti equations 

A2(x; t) ± A'(x; r) = 2k'^sn'^{x ± r) - /c^ + £(,-) . (3.9) 

Indeed, from factorizations p.7p it follows that the T){x;t) and P^(x;r) intertwine the Hamil- 
tonians H{x^) and H(x-), 

V{x;t)H{x+) = H{x^)V{x;t), V\x;t)H{x-) = H{x+)V\x;t) , (3.10) 

and interchange the eigenstates of the superpartner systems, 

P(x;r)M/^(x+)= J-|'(a,r)M/^(x_), V\x;T)nix.) = -T^{a, -t)^'^{x+) . (3.11) 

The second relation in (13. lip follows from the first one via a substitution r — t- — r. A complex 
amplitude, I'^{a,T) = e=^*'^^("'^)A^^(a, r), is given by 

^■D, , r / . X TT \ 1 e(2r±a)e(o) 

It satisfies (7"^(a,r))* = J-'^(q,t) = — -F^(a, — r). Its modulus may be presented in a form 
A4'^{a,T) = E{a) + e(T) , where £'(a) for the valence and conduction bands is given by Eqs. 
(j2.5p and (j2.6p . This agrees with Eq. (j3.7p . Notice that the modulus is an even in r function, 
A4^{a,T) = Ai^{a, — r), which is nonzero except for the lower edge states of the valence band 
(E = 0) in the case t = + n)K. A phase is well defined for ^ 0, and satisfies a relation 

^i<fi'^{a-T) _ _g-j'y3^ (a.r) ^ (3.13) 
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E=l 




Figure 2: The self-isospectral potentials V± = 2A;^sn(x-i-) — A;^ are shown together with the 
edges of the valence (0 < < A;'^) and conduction {1 < E < co) bands. V± have maxima at 
X = =Fr + (2n + 1)K and minima at x = =ft + 2nK. Here A;^ = 0.75, K = 2.16, and r = 0.8. 



It can be presented in a form 



gi<P^(«,r) ^ -sign (ns 2r) exp -2i (^^(a) - ^) r + i(^e(a, t)J , (3.14) 

where sign(.) is a sign function, and (^0(a,r) is a phase of 0(2r + a), 999(0, r) = 
Im Z(u) du) , see Eq. (ICT . Particularly, for the edge states (i = 1,2,3), Eq. (IXT2]1 

gives P(x;r)V'i(x+) = (r)V'i(x_), T-)^.(a;_) = jr^(r)V'i(x+), where 

■0i(x)=dnx, cnx, snx, J^f (r) = — cn2rns2r , — dn2rns2r, — ns2r, (3.15) 

and so, 

Mf{T) = ^J^, y^k'^ + e{T), V^ + ^ir), (3.16) 

and e^^i^^'^^ = — sign (cn 2r ns 2t), — sign(ns2r), — sign(ns2r). 

As a consequence of intertwining relations (j3.10p . first order matrix operators 

are the integrals of motion for system (j3.8p . Integrals (j3.17p correspond here (up to a unitary 
transformation of sigma matrices) to the first order operators Ga in section [TJ Operator T = cr^ 
is a trivial integral for (j3.8p . [E,"?^] = 0, that anticommutes with Sa, a = 1,2, {T,Sa} = 0, and 
classifies them as supercharges. Bosonic, Ti, and fermionic, Sa, operators satisfy then the N = 2 
supersymmetry algebra, 

{Sa,Si,] = 26ab{n + e{T)), [n,Sa] = Q. (3.18) 
In correspondence with (13. lip and ()3.13p , the eigenstates of the supercharge Si are 
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Since e(T) > for r 7^ + n)K, n G Z, the first-order super symmetry ()3.18p R is dynamically 
broken in general case. It is unbroken however for r = (n + i)K by virtue of e{{^ + n)K) = 
0. For these values of the shift parameter, the supercharges Sa annihilate the ground states 
dn(x + (n + ^)K) and dn(x — (n + ^)K) of the super-partner systems H{x + (n + ^)K) and 
H{x — {n + ^)'K). Notice that with variation of the shift parameter r 7^ nK, which simultaneously 
governs the scale of the supersymmetry breaking £{t), the spectrum of the second order system 
(j3.8p does not change. Each of its two super-partners has the same spectrum as a non-shifted 
Lame system ()2.ip does. Therefore, each energy level inside the valence, < E < k'^, and 
conduction, 1 < E < 00, bands is fourth- fold degenerate in accordance with the existence of the 
two Bloch states, and of the form (j2.3p for each subsystem, see Eq. (I3.19p . 

We have a two-fold degeneration at the edges E = 0, E = k''^ and = 1 of the valence 
and conduction bands in the spectrum of the super symmetric system Ti. Bosonic, and 
fermionic, states are defined as eigenstates of the grading operator T = cjs, E^'^^-' = ib^'^^\ 

and have the general form ^'(+) = (\I'(a;-|-), 0)^ and ^f^") = (0, where T means a 

transposition. In summary, we see that in both the broken and unbroken cases, the Witten index, 
which characterizes the difference between the number of bosonic and fermionic zero modes, is 
the same and equals zero. 

Eor r 7^ + n)K (when e(r) 7^ 0) super symmetric relations (|3.18p look differently from a 
usual form of superalgebra in supersymmetric quantum mechanics. A simple redefinition of the 
matrix Hamiltonian ()3.8p . % ^ % = % + £(''"); will correct the form of super algebraic relations, 
but will not change the conclusions on a broken (for r 7^ (^ + n)K) form of the supersymmetric 
structure that we have analyzed. We shall return to this point in the discussion of the peculiar 
supersymmetry of the first order Bogoliubov-de Gennes system in section [71 

The described degeneracy of the energy levels in both, broken and unbroken, cases is unusual 
for N = 2 supersymmetry. We will show that additional nontrivial integrals of motion may be 
associated with this peculiarity of the self-isospectral supersymmetric system (13. Sp . To identify 
such integrals, in the next section we investigate the function A(x;r) in greater detail. 



4 Superpotential 

Being the logarithmic derivative of F{x;t), see Eq. ()3.6p . the superpotential A(x;r) may be 
written with the help of (IB.lip . (IB.14P in terms of Jacobi's Z, or and H functions, 

A(.;r) = z(r) + Z(.-) - Z(.,) = l|ln ( g., "'g,,^, ) ■ M 

The addition formula (IB.6P for the Z function gives another, equivalent representation 

A(x;r) = ?(r) + A;^sn2rsn (x_) sn (x+) . (4.2) 

Functions z(r) and ?(r) are defined in (j3.2p . (j3.3p . Yet another useful representation for the 
superpotential may be derived from (j4.2p . 



snx_cn2;_dnx_ +snx+cnx+dnx+ 
sn^X-i- — sn^a;_ 

Having in mind relations ()3.10p . ()3.7p and (13. 9p . in what follows we treat x as a variable 
and r as a shift parameter. A(x; r) is an elliptic function in both its arguments with the same 



This refers to the order of the polynomial in H that appears in the anticommutator of the supercharges. 
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periods 2K and 2zK'. It is an even in x and odd in the r function with respect to the points 
0, K (modulo periods), A(-x; r) = A{x; r), A(K -x;t) = A(K + x; r), A{x; -r) = -A(x; r), 
A(x; K-r) = -A(x; K+r). It also obeys a relation A(x+K; r+K) = A(x-K; r+K) = A(x; r). 
In r = 0, K the function undergoes infinite jumps. 

Being the elliptic function in x, A(x;t) obeys a nonlinear differential equation 

A'^ = A^ + 252(r)A2 + 5i(t)A + 5o{t) , (4.4) 

where (52 (r) = 1 + k"^ - 3ns'^2T, <5i(r) = Sns^ 2r cn2r dn2r, and So{t) = -3ns'^2r + 2(1 + 
A;^)ns^ 2r + k''^. As a consequence of (j4.4p . it also satisfies the nonlinear higher order differential 
equations 

A" = 2 A3 + 2^2 (r) A + hiir) , A'" = 2 A' (SA^ + 62{t)) . (4.5) 

Making use of (|4.ip . one finds the relation 

A(x + r + A; A) - A(x + A; r + A) + A(x; r) = g{T, A) . (4.6) 

Function g{T,X) = (^{t) + <;^(A) — <j(t + A) + A;^sn2r sn2Asn2(r + A) has symmetry properties 
g{T, A) = g{X, t) = g{T, —A — r) = —g{—T, —A), and may be written as 

g{T, A) = ns 2r ns 2A ns 2(t + A) [1 - cn 2r cn 2A cn 2(t + A)] . (4.7) 

For a particular case A = K/2, to be important for non-periodic limit, 

g {t,^K) =C{t) , C(t) = ns2rnc2rdn2T. (4.8) 

Notice that g{T, A) takes nonzero values for all real values of its arguments El- Equation (14. 6p is a 
kind of addition formula for elliptic function A(x;r). Differentiating (|4.6p in x and using Ricatti 
equations ()3.9p . we obtain a relation 

A'(x + r + A; A) - A(x + A; r + A)A(x + r + A; A) = 
(A2(x; t) + A'(x; r) + 52(r)) - ^(r, A)A(x; r) + G(r, A) , (4.9) 

where G(r, A) = ^ [l + P + c/2(r, A) - ns^ 2r - ns^ 2A - ns^ 2(t + A)] = 0. 

In conclusion of this section we note that functions 5a{T), a = 0,1, 2, can be given a physical 
sense by expressing them in terms of the band edges energies and of e(r) : S2{t) = — (£'^-|-£^|-|-£'|), 
6iiT) = -2^, 5o{t) = -52{t) - 2(^1^2 + ^1^3 + ^2^3), where E,{t) = E., + e(r), Ei = 0, 
E2 = fc'^ and £"3 = 1. Particularly, 5i measures a velocity with which a scale of supersymmetry 
breaking changes as a function of the shift parameter. Notice also that the first equation in (j4.5p 
has a form of a modified Ginzburg-Landau equation, see 03], which corresponds here to a gap 
equation for the real condensate field in the kink-antikink crystalline phase in the Gross-Neveu 
model with a bare mass term, see PE]. At r = (i +n)K we have 5i = 0, and superpotential 
A(x) satisfies the nonlinear Schrodinger equation, the lowest nontrivial member of the modified 
Korteweg-de Vries hierarchy [44]. This homogenisation of the second order nonlinear differential 
equation can be associated with restoration of the discrete chiral symmetry in (jl.2p at mo = 0. 

^It takes zero values at some complex values of the arguments, for instance, C(iK ± f K') — 0. 
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5 Higher order integrals and nonlinear superalgebra 



Now we are in a position to identify higher order local intertwining operators and integrals of 
motion for the system %. First, we find the second order intertwining operators. Changing 
r — )• —A and shifting the argument x— T-x + r + Ain the first relation from (I3.10p . we obtain 

V{x + T + A; -\)H{x + t) = H{x + t + 2\)V{x + r + A; -A). (5.1) 

Multiplying (15. ip by V{x + A; r + A) from the left, and using once again (13.10p on the right hand 
side, we obtain an intertwining relation 

B{x] T, X)H{x+) = H{x-)B{x] T, A) . (5.2) 

It is generated by the second order differential operator 

Bix; T, A) = V{x + A; r + X)V^{x + r + A; A) , (5.3) 

which is defined for A , r + A ^ nK. For adjoint operator we have B^{x; r, X)H{x — r) = H{x + 
t)B^ {x;t, X). In accordance with ()5.ip . the second order intertwining operator (j5.3p shifts the 
Hamiltonian's argument first for 2A and then for — 2(t + A). Equivalent representation of the 
operator (15. 3p is 

B{x; T, A) = -yix; r) - g{T, X)V{x; r) , (5.4) 

3^(x;r) = ^-A(x;T)^-i(A2(x;r) + A'(x;r)+52(r)) , y\x- r) = y {x; -r) . (5.5) 

We have used here Eq. (j4.6p . So, the dependence of B{x;t,X) on A is localized only in the 
x-independent multiplier g{T,X), see Eq. (j4.7p . 

From Eqs. (15. 3p and (I3.10p it follows that at r = the second order intertwining operators 
B{x;t,X) and B^x;t,X) reduce, up to an additive term e(A), to the isospectral superpartner 
Hamiltonians, B{x; 0, A) = H{x) + e(A), | B^^ix; 0, A) = H{x + 2A) + e{X). 

Forgetting for the moment on the r = case, from the viewpoint of intertwining relation 
(j5.2p . one could conclude that the parameter A has a "gauge- like" , non-observable nature. Such 
a conclusion, however, is not correct. We will return to this point later. 

Since g{T, A) is nonzero for real r and A, operator y{x; r), unlike B{x; r. A), is not factorizable 
in terms of our first order intertwining operators (with real shift parameters) Nevertheless, it is 
the second order intertwining operator as well as B{x; r, A). It can be presented as a linear combi- 
nation of the second and first order intertwining operators, y{x; r) = —B{x; r, A) —g{T, X)T>{x; r), 
and also may be used together with the first order operator P(x; r) to characterize the system. 
At the end of this section we shall discuss the peculiarities associated with such an alternative. 

Having in mind a non-periodic limit we discuss later, it is convenient to fix A = K/2, and 
introduce a notation A{x;t) = B (x;r, ^K), i.e. 

A{x; t) =V{x + ^K;t + ^K) V^x + t + ^K; ^K) = -y{x; r) - C{t)V{x; t) , (5.6) 

where C(r) is defined in Eq. (j4.8p . Employing the properties of y{x;T) and 'D{x;t) under 
hermitian conjugation, from ()5.6p one finds A^{x;t) = A{x;—t), and then a representation 

®One could conclude that Eq. H5.4|) contradicts to this relation since g{T,\) diverges at r = 0, and operators 
V{x;t) and 3^(x; r) are not defined for r — 0. Eq. (|5.4p correctly reproduces this relation by treating r = as a 
limit r — > 0, and employing addition formulae HA.6|I for Jacobi elliptic functions. 

''It can be factorized in terms of our first order Darboux operators D in special cases of r = (i + n)K. Such a 
factorization corresponds to complex values of the shift parameters, see a discussion below in this section. 
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alternative to ([52]) is obtained, A{x; r) = P(x - r + ^K; ^K)'D'f{x + ^K; -r + ^K). Unlike 
the operators V{x;t) and y{x;T), the ^(x;r) is well defined at r = and reduces just to a 
non-shifted Hamiltonian, ^(x;0) = J^{x;Q) = H{x). Notice, however, that unlike P(a;;r), it is 
not defined for t = + n)K. 

Second order intertwining operator of the most general form (j5.3p may be presented in terms 
of the intertwining operators A{x] r) and V{x] r), B{x; r, A) = A{x] r) + (C(t) — g{T, A))P(x; r). 

Because of Eq. (j5.2p . the self-isospectral system possesses (for r 7^ (| + n)K) the second 
order integrals 

"^'T')- * = 

to be nontrivial for r 7^ nK and independent from the first order integrals p.l7p . 
With some algebraic manipulations, we find 

A\x]t)A{x]t) = H{x+)[H{x+) + q{t)] , where £.(r) = fc'^sn^ 2r nc^ 2r . (5.8) 

A similar relation is obtained from (j5.8p by a simple change r — )• — r, ^(x; r)^^(x; r) = 
H{x-) [H{x-) + i?(T)] . cf. relations in (j3.7p for the first order intertwining operators. 

The intertwining second order operator A{x] r) annihilates the lower energy state dn (x + r) 
of the system H(x + t). Another state annihilated by it is 

fix, r) = dn (X + r) f + ^.^'^ + rf^ , (5.9) 
J dn (ti + r) 

and we have f{x + 2K,t) = exp [2Kz(r + ^K)] f{x,T). Function (15. 9p for r / is unbounded 
and describes therefore a non-physical eigenstate of H[x + t) from the lower forbidden band with 
energy E = —q{t) < 0, see Eq. (15. Sp . At r = 0, function (|5.9p reduces to E(x + K)dnx that 
corresponds to a nonphysical state of H{x) of zero eigenvalue. 

Like the first order operator D^x; r), A{x; r) transforms the eigenstates of H{x + t) into those 
ofH{x-T), 

A{x;t)^1{x+) = J-^(a,r) ^^(x_) , (5.10) 

where 

T^{a, t) = e^'^'^^'^'^^M^ia, r) , M-^{a, r) = y^E{a){E{a) + g{T)) . (5.11) 

The modulus and the phase of the complex amplitude J-^{a,T) are expressed in terms of those 
for the first order intertwining operator by employing Eqs. (j5.ip . (j5.6p and (|3.1ip . 

M-^{a, r) = {a, r + ^K) (a, ^K) , ip-^{a, r) = (^^(a, r + ^K) - V9^(a, ^K) . (5.12) 

A phase (f^{a,T) G M has, unlike (j3.13p . a property e*'^'^^"'"'^) = e"*'^'^^"''^) due to a relation 
^^(x;r) = ^(x; — r) to be different in sign from that for the first order intertwining opera- 
tor, V^{x;t) = —V{x;—t). For the edge band states, particularly, we have A{x^,T)^pi{x^) = 
JV^(r)'0j(x_), A^{x■,T)^l'i{x^) = Ff{T)il)i{x+), where J^fir) = 0, A;'^nc2r, dn2rnc2r, i = 
1, 2, 3, cf. (|3.15p . The eigenstates of the integral Qi, see (|5.7p . have a form similar to that for Si, 

Two relations are valid for the first and second order intertwining operators, 

V\x; t)A{x; t) = V{x+) - C{t)H{x+) , V{x; t)A\x; t) = -P(x_) - C{t)H{x^) . (5.14) 
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Here V{x^) = V{x + t) is an anti-hermitian third order differential operator 

d 

= -— ^ + (l + A:^ — 3fe^sn^x+) 3/i;^snx+ cna;+ dnx+ . (5.15) 

dx-^ ' dx 

Notice tliat like the Lame Hamiltonian, the operator (j5.15p is well defined for any value of the 
shift parameter r. Two related equalities may be obtained from (j5.14p by hermitian conjugation. 
Making use of intertwining relations (|3.10p . (j5.2p . we find that li{x + r) commutes with 
r)^(j;; r), and, therefore, Vix + r) is an integral for the subsystem li[x + r). For self- 
isospectral supersymmetric system 7^ we have then two further, third order hermitian integrals 



Li = -idiag L2 = usLi . (5.16) 

Operator P(x) is a Lax operator for the periodic one-gap Lame system H{x), see [38| [39]. 
The following relations that involve the operator V{x^) are valid, 

P(x; t)V{x + t) = A{x; t) [H{x+) + e(r)] + C(r)P(x; t)H{x+) , (5.17) 

A{x; t)V{x+) = -V{x; t)H{x+) [H{x+) + q{t)] - C{t)A{x; t)H{x+) , (5.18) 
-V^{x+) = P{H{x+)), P{H) = H{H -k'^){H -I). (5.19) 



The third order polynomial P{H) is the same spectral polynomial of the Lame system that arose 
before in (j2.7p and in differential dispersion relation (j2.1ip : it turns into zero when acts on the 
edge states with energies Ei = 0, A;'^, 1. Since the third order differential operator 7^(x+) is an 
integral of motion for i7(x+), relation ()5.19p means that the edge states dnx+, cna;_|. and sna;_|_ 
form its kernel [39]. The spectral polynomial is a semi-positive definite operator, while V{x) is 
an anti-hermitian operator. Its action on physical Bloch states (j2.3p should reduce therefore to 
■iii-^P{E{a)). The phase cannot change abruptly within the allowed bands. To fix correctly 
the sign, one can consider a limit — )• 0, in which Lame system (|2.ip reduces to a free particle, 
an integral V{x) reduces to a third order operator d^ /dx'^ + d/dx, forbidden zone k''^ < E < 1 
disappears, Bloch states transform into the plane wave states, whereas the edge states dnx, cnx 
and snx reduce, respectively, to 1, cosx and sinx with energies E = 0, 1 and 1. Summarizing 
all this, one finds that the operator (|5.15p acts on the physical Bloch states (j2.3p as follows, 

Vix)^1{x) = TME)y/p{E{a))^1{x) , (5.20) 

where, as in (j2.7p and (j2.1ip . ri{E) = —1 for valence and +1 for conduction bands Relation 
()5.20p means, particularly, that the Lax operator is not reduced just to a square root from the 
spectral polynomial since Hamiltonian does not distinguish index it. This is a true, nontrivial 
integral of motion that is related with the Hamiltonian H by polynomial equation (I5.19|) Eq. 
()5.19p corresponds to a non-degenerate spectral elliptic curve of genus one associated with a 
one-gap periodic Lame system [35] . 

Applying the first relation from (I5.14p to a physical Bloch state (a;+) and using an equality E{E ~\~ q{t)){E + 
e(T)) — P{E) + C {r)E , we obtain the Pythagorean relation for a rectangular triangle with legs C(r)iJ(a) and 
^P{E{a)), V-P(-B(a))+ C2(T)S2(a)e*(^''("'"+f » ^ iriy' P{E{a)) + C{T)E{a). 

^This corresponds to Burchnall-Chaundy theorem [45] that underlies the theory of nonlinear integrable systems 
[35| . It asserts that if two ordinary differential in x operators A and B of mutually prime orders I and m do 
commute, they obey a relation P{A,B) — 0, where P is a polynomial of order m in A, and of order I in B. 
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Let us discuss now the superalgebra generated by the zero, CT3, first, Sa, second, Qa, and third. 
La, order integrals of motion of the self-isospectral system Ti. The operator T = commutes 
with La and anti-commutes with Qa, and so, classifies them, respectively, as bosonic and fermionic 
operators. Using the displayed relations for the operators V, A and V as well as those obtained 
from them by a hermitian conjugation and by a change t ^ —t, Eq. (|3.18p is extended by the 
ant i- commutation relations of the integrals Sa with Qa, and the commutation relations of Sa and 
Qa with La- We arrive as a result at the following superalgebra for the self-isospectral system 



(j3.8p with Z2 grading operator F = 1T3, 

{Sa,Sa} = 25ab{n + eiT)), {Qa, Qb} = 25abn{n + g{T)) , (5.21) 

{Sa, Qb} = 2 {-6ab C{T)n + eab ^l) , (5.22) 

[Li, Sa] = [Li, Qa] = [Li,L2] = , [L2, Sa] = 2i{SaC{T)n + Qa{n + e(T))) , (5.23) 

[L2, Qa] = -2i{San{n + q{t)) + QaC{T)V) , (5.24) 

ksjSa] = -2ieabSb, [0-3, Qa] = -2ieabQb, [c^3,La] = 0, (5.25) 

[n, as] = [n, Sa] = [n, Qa] = [n, = o . (5.26) 



We have here a nonlinear superalgebra, in which Li (that is a Lax operator for H) plays a role 
of the bosonic central charge, and is treated as one of its even generators in correspondence 
with Z2 grading relations [as, 0-3] = [crs,^] = [as. La] = and {a3,Sa} = {cr3,Qa} = 0. 

Since Li commutes with Sa and Qa, the eigenstates (j3.19p and (|5.13p of and Qi are 
simultaneously the eigenstates of Li, 

where A = S*! or Qi, rj is the same as in ()2.1ip and (I5.20p . and P{a) = P{E{a)). Note that 
unlike Si and Qi, Li distinguishes the index it. 

In correspondence with the discussion related to (|5.9p . the Qa, a = 1,2, annihilate the two 
ground states of zero energy, dn (x + r) and dn (x — r) , while other two states from their kernel 
are non-physical. These supercharges are not defined, however, for r = (^ + n)K, which are the 
only values of the shift parameter when the N = 2 supersymmetry associated with the first order 
supercharges Sa is not broken. Therefore, when the first and the second order supercharges are 
simultaneously defined [for r / (| + n)K, nK], the supersymmetry generated together by Sa and 
Qa is partially broken. 

One could construct, instead, the second order supercharges, Q^, on the basis of the inter- 
twining operators 3^(x; r) and y'^{x;T). According to (|5.6p . they are related to Qa as 

Q^a = -Qa - C{T)Sa ■ (5.28) 

The corresponding super-algebra with Qa substituted for Q^ will have then a form similar to 
that we have discussed, with the change of some of the corresponding (anti)-commutators for 

{Ql, Q^} = 26ab{nin + gir) - C\t)) + £(r)C2(r)) , (5.29) 

{Sa,Q^} = -2{dabasC{T)£{T)+eabLi) , (5.30) 

[L2, Sa] = -2i{SaCiT)e{T) + Q^iH + e(r))) , (5.31) 
[L2, Q^] = 2i{San{n + gir) + e{T)C(r) - C\r)) + Q^e(r)C(T)) . (5.32) 
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The second order supercharges Q-^, like Sa, are well defined at r = + n)K but not defined 
for T = nK. Analyzing the roots of the polynomial in the right hand side of (j5.29p . one finds 
that the kernels of Q^, a = 1,2, for r / (| + n)K are formed by non-physical states. In the 
exceptional case t = + n)K, for which the supercharges Qa are not defined, the polynomial in 
(j5.29p reduces to the second order polynomial 

PQy{'H) = {n-k'^){n-i). (5.33) 

In correspondence with this, the zero modes of the operators 3^(2;; ^K) and y^{x; ^K) = 
y{x; — ^K) are, respectively, the physical edge states cn {x + ^K), sn (x + ^K) and cn (x — ^K), 
sn (x — ^K). This property reflects a peculiarity of the case r = (| + n)K in another aspect. 
In accordance with footnote 5, function g{T, A) in (j5.4p turns into zero at A = i(K + iK'). The 
second order operator y{x; ^K) factorizes then either as y{x; |K) = —V{x + ^(K + iK'); K + 
\iK')V\x + K + iiK'; i(K + iK')), or in alternative form obtained by the change of i for —i. 
These two factorizations can be presented equivalently as 

y{x- iK) = (ns (x - iK)^sn (x - iR)) (cn (x + iK)^nc (x + iR)) , (5.34) 
y{x- iK) = (nc (x - iK)|:cn (x - ^K)) (sn (x + iK)^ns (x + ^K)) . (5.35) 

From here we see that the particular case of the half period shift of the super-partner systems 
is indeed exceptional. In this case not only the N = 2 supersymmetry associated with the first 
order supercharges Sa is unbroken (when zero modes of Sa are the ground states that form a 
zero energy doublet), but all the other edge states of the energy doublets with E = k'"^ and 
E = 1 correspond to zero modes of the second order supercharges Q^. Then the third order 
spectral polynomial -P(^) = 'H{'H — k''^){T-L — 1) is just a product of the first and the second 
order polynomials which correspond to the squares of the first, Sa, and the second, Q^, order 
supercharges. In this special case the (anti-)commutation relations (I5.30p . (15.311) . (|5.32p also 
simplify their form, {Sa,Q^} = -2e,feLi, [L2,Sa] = -2iQyn, [L2,Q^] = 2iSaPQy{n). We also 
have 

SaQ^ = -Q^Sa = -iL2 , SaQ^ = Sa = "Li , (5.36) 

where there is no summation in index a, and b ^ a. This is in conformity with the above men- 
tioned factorization of the spectral polynomial. However, since Qa does not annihilate the ground 
states dn (x -|- ^K) and dn (x — ^K) (which are transformed mutually by the intertwining opera- 
tors y{x; ^K) and y^{x; ^K)), we conclude that nonlinear supersymmetry of the self-isospectral 
system also is partially broken at r = (^ + n)K0. 

In the next section we will see that another peculiarity of our self-isospectral system is that 
the choice F = (T3 is not unique for identification of the Z2 grading operator : it also admits 
other choices for F, which lead to different identifications of the integrals (T3, Sa, Qa and La as 
bosonic and fermionic operators. This results in the alternative forms for the superalgebra. Each 
of such alternative forms of the superalgebra makes, particularly, a nontrivial relation (j5.19p to 
be 'visible' explicitly just in its structure, unlike the case with F = C73 that we have discussed. 
We also will identify the integrals of motion which detect the phases in the structure of the 
eigenstates of the operators Sa and Qa- 

Cf. this picture as well as that for r 7^ (i +n)K, which we discussed above, with the picture of supersymmetry 
breaking in the systems with topologically nontrivial Bogomolny-Prasad-Sommerfield states !46] . 
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6 Nonlocal Z2 grading operators 



Let us introduce the operators of reflection in x and r, TZxTZ = —x, TZtTZ = r, TZ^ = 1, TtT = — r, 
TxT = X, = 1. They intertwine the superpartner Hamiltonians, TZH^xj^) = H{x )lZ, 
TH{x+) = H{x-)T, and we find that the self-isospectral supersymmetric system (I3.8P possesses 
the hermitian integrals of motion 

TZai, Tdi, TZa2, Ta2, TlTa^, TZT . (6.1) 

Like for (T3, a square of each of them equals 1. From relations 

TZ'D{x; t) = V\x]t)TZ, TZA{x-t) = A\x;t)1Z, llV{x+) = -V{x-)ll , (6.2) 
TV{x-t) = -V\x;t)T , TA{x;t) = A\x;t)T , TV{x+) = V{x^)T , (6.3) 

it follows that IZ and T intertwine also the operators of the same order within the pairs iT){x] r), 
T>'^{x;t)), {A{x;t), A'^{x;t)), and (7^(x+), V{x-)). As a result, each of the nonlocal in x or r, 
or in both of them, integrals of motion (|6.ip either commutes or anti-commutes with each of the 
nontrivial local integrals Sa, Qa and La- Then each integral from (16. ip also may be chosen as 
the Z2 grading operator for the self-isospectral system (jS.Sp . Corresponding Z2 parities together 
with those prescribed by a local integral (T3 are shown in Table [H Z2 parities of the second 
order integrals Q^, defined in ()5.28p . are also displayed; the equality C(— r) = — C(t) has to be 
employed in their computation. Notice that Q^, a = 1,2, always has the same Z2 parity as the 
Qa with the same value of the index a. 



Table 1: Z2 parities of the local integrals. 
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A positive Z2 parity is assigned for the Hamiltonian H by any of the integrals (j6.ip . Then 
for any choice of the grading operator presented in Table [H four of the eight local integrals , 
T-L, Sa, La and Qa or Q^ are identified as bosonic generators, and four are identified as fermionic 
generators of the corresponding nonlinear superalgebra. The superalgebra may be found for each 
choice of F from the set of integrals (16. ip by employing the quadratic products of the operators 
V, A and V that have been discussed in the previous section. Alternatively, some of the (anti)- 
commutators may be obtained with the help of the already known (anti)-commutation relations 
and relations between the generators that involve a^. For instance, [5i,Qi] = i<Jz{S\,Q2}- As 
an example, we display the explicit form of the superalgebraic relations for the choice F = TZT, 

{Sa,S^} = 25ab{n + e{T)), {Sa,Li} = 2eai>{Qb{n + e{T))+C{T)Si,n), {Sa,L2} = Q, (6.4) 
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{Li, Li} = {L2, L2} = 2P{n), {Li, L2} = 2a3P(7^), (6.5) 
[Qa,Sb] = 2i{-6abL2 + eabCiT)a3n), [Qi,Q2] = -2iam{U + q{t)), (6.6) 

[Q„,Li]=0, [Qa,L2]=2i[C{T)Qan + San{n + Q{T))), (6.7) 

which should be supphed by the commutation relations (|5.25p and (j5.26p . PiU) in (j6.5p is the 
spectral polynomial, see (j5.19p . 

A fundamental polynomial relation (j5.19p between the Lax operator and the Hamiltonian, 
that underlies a very special, finite-gap nature of Lame system ^1, does not show up in the 
superalgebraic structure for a usual choice of the diagonal matrix as the grading operator F, 
but is involved explicitly in the superalgebra in the form of the anticommutator of one or both 
generators Lq, a = 1,2, when any of six non-local integrals (16. ip is identified as F. 

Note that for F = TZT as well as for any other choice of the grading operator that involves the 
operator 7~, the constant C(r) anticommutes with the grading operator and should be treated as 
an odd generator of the superalgebra. As a result, the right hand side in the second anticommu- 
tator in (j6.4p is an even operator, while the right hand side in the first (second) commutator in 
(j6.6p (in (j6.7p ) is an odd operator as it should be. 

By employing Eq. (15.28p . one can rewrite superalgebraic relations (16. 4p . (16. 6p and (|6.7p in 
terms of the integrals Q^, which, unlike Qa-, are defined for r = (^ + n)K. We do not display 
them here, but write down only a commutation relation 

[Sa, Q^] = 2i{6abL2 + a3eabC{T)e{T)) , (6.8) 

which we will need below. The form of such a superalgebra simplifies significantly at r = (^ + 
n)K in correspondence with a special nature that the integrals Sa and acquire in the case. 
Particularly, one finds 

{Sa, Sb} = 26ab'H , {Sa, L^] = -2eabQ^n , (6.9) 
[Qy,Sb]=2i5abL2, [Qy,Qy] = -2iasPQy{n), [L2,Qy] = 2iSaPQy{'H) . (6.10) 



All the integrals (16. ip including as but excluding TZT may be related between themselves 
by unitary transformations, whose generators are constructed in terms of the grading operators 
themselves. For instance, Ua^W = IZai = ffs, U = U"^ = U^^ = -^{(^3 + T^ci). Being con- 
structed from the integrals of motion, such a transformation does not change the super symmetric 
Hamiltonian Ti. On the other hand, if we apply it to any nontrivial integral, the transformed 
operator still will be an integral. Particularly, its application to the integrals Si and Qi gives 

S = ^7^c^25'l = diag (nV{x; r), -nv\x; r)^ , Q = ina2Qi = diag (nA{x; r), -TZA^x; t)J . 

(6.11) 

These are nontrivial hermitian nonlocal integrals of motion for the self-isospectral system ()3.8p o. 
Eq. (16. lip has a sense of Foldy-Wouthuysen transformation that diagonalizes the supercharges 
5i and Qi. The price we pay for this is a non-locality of the transformed operators. 

Multiplication of (16. lip by the grading operators gives further nonlocal integrals, particularly, 
a^S and a^Q. Since the both operators (|6.1ip are diagonal, the Lame subsystem H{x^) may be 
characterized, in addition to the Lax integral V{x^), by two nontrivial nonlocal integrals 

S = nV{x;T), Q = nA{x;T). (6.12) 

^^In a generic situation the spectrum of a one-dimensional periodic system has infinitely many gaps [35j . 
^^Notice that the (f +f )-dimensional GN model has a system of infinitely many (nonlocal) conservation laws. 
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In correspondence with relations P^(x;r) = — — r) and A^{x;t) = A{x;—t), another sub- 
system H{x-) is characterized then by the integrals of the same form but with r changed for 
— T. The operator T = TZT is an integral for the subsystem H{x-^-) [as well as for subsystem 
H{x^)]. It can be identified as a Z2 grading operator that assigns definite Z2 parities for non- 
trivial integrals of the subsystem H{xj^). Namely, in correspondence with (j6.2p and (j6.3p . the 
integrals — i7^(x+) and S are fermionic operators with respect to such a grading, while Q should 
be treated as a bosonic operator. Multiplying fermionic integrals by iV and bosonic integral by 
r, we obtain three more integrals for H{xj^). It is not difficult to calculate the corresponding 
superalgebra generated by these integrals. Let us note only that since the described supersym- 
metry may be revealed in the subsystem H{x^) (or, in H(x-)), it may be treated as a bosonized 
super symmetry, see [171 ETJ [38] . 

Let us return to the question of degeneration in our self-isospectral system. This will al- 
low us to observe some other interesting properties related to the nonlocal integrals (|6.ip . 
Let us take a pair of mutually commuting integrals Si and Li. They can be simultaneously 
diagonalized, and for their common eigenstates we have •S'l^'^ ^ = eAi^{a,T)"^'^ ^ and 

^i^±,Si,e = =Ff?(a)\/-P(a) ^±,Si,e^ see Eqs. ([3T9]) and ([OT]) . We can distinguish all the four 
states by these relations for any value of the energy within the valence and conduction bands, 
and each two doublet states for the edges E = 0, k'"^, 1 of the allowed bands when r 7^ -|- n)K. 
However, in the case of r = (| -|- n)K, the two ground states of zero energy are annihilated by 
the both operators Si and Li, and cannot be distinguished by them. In this special case the 
operator as commutes with Si and Li on the subspace E = 0, and may be used to distinguish 
the two ground states. It is necessary to remember, however, that as does not commute with Si 
on the subspaces of nonzero energy. 

There is yet another possibility. According to Table [TJ the local integrals Si and Li commute 
with the nonlocal integral To"2- We find then 

7"^2^^,5i,e = ^ee^^^"^"'^^^±,5i,., (6.13) 

where we used relation (j3.14p . The operator Ta2 detects therefore the phase in the structure of 
the eigenstates of 5*1. By comparing two supersymmetric systems with the shift parameters r 
and T -|- K, and by taking into account the 2K-periodicity of the function in (j3.12p and the 
2K-anti-periodicity of snit, we get from ([TOD that e<¥'^(".^+K)-v'^(a,r)) ^ ^±^{a)T ^ ^^^^^ ^^le 
integral 'Ta2 makes, particularly, the same job as a translation for the period operator (which is 
also a nonlocal integral for the system) : it allows us to determine an energy-dependent quasi- 
momentum. Finally, in the case of zero energy {a = K. + iK'), treating r = (^-|-n)K as a 
limit case, one can also distinguish two ground states in the supersymmetric doublet by means 
of (I6T3]1 . 

Instead of 5i, Li and Ta2, we could choose the triplet ^2, Li and Tai of mutually commuting 
integrals, see Table [H The states within the supermultiplets can be distinguished also by choosing 
the triplets of mutually commuting integrals (Qi, Li, Tai), or {Q2, Li, T1T2). For the two latter 
cases, the doublet of the ground states is annihilated by Qa and Li for any value of the shift 
parameter r (excluding the case t = + n)K when Qa are not defined), but the corresponding 
integrals Tai or Ta2 do here the necessary job of distinguishing the states as well. 

The integrals TZai and IZTa^ act on the eigenstates of Si, with which they also commute, 
as lZai'^1s^^^{x,T) = -ee=^*'^'^("'^)^'^^ ^^^^(x, r), 7^r^T3*^ <j^^^(x, r) = -'^^ g^ J^x,T). These op- 
erators interchange the states with -|- and — indexes, and anti-commute with the integral Li. 
The edge states, which do not carry such an index, are annihilated by Li, so that there is no 
contradiction with the information presented in Table [TJ 
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In conclusion of this section we note that the Witten index computed with the grading oper- 
ator identified with any of the six nonlocal integrals (I6.ip is the same as for a choice F = (T3, i. e. 
/^w = 0. 

7 Super symmetry of the associated periodic BdG system 

Till the moment we have discussed the self-isospectrality of the one-gap Lame system with the 
second order Hamiltonian. Though we have shown that its supersymmetric structure is much 
more rich than a usual one, from the viewpoint of the physics of the GN model it is more natural 
to look at the revealed picture from another perspective. 

Let us take one of the first order integrals Sa of the self-isospectral Lame system, say Si, and 
consider it as a first order, Dirac Hamiltonian. In such a way we obtain an intimately related, 
but different physical system. Unlike the second order operator the spectrum (I3.19P of Si 
depends on r. We get a periodic Bogoliubov-de Gennes system with Hamiltonian HsdC = 'S'l- 
The interpretation of the function A(j;; r) changes in this case : this is a Dirac scalar potential in 
correspondence with a discussion from section [TJ In dependence on a physical context, it takes a 
sense of an order parameter, a condensate, or a gap function. 

The r-dependent spectrum of such a BdG system consists of four or three allowed bands 
located symmetrically with respect to the level £ = see Figure [3l Interpretation of the bands 
also changes and depends on the physical context. For r 7^ + n)K, the positive and nega- 
tive 'internal' bands are separated by a nonzero gap A£^(r) = 2-^/ e(r) = 2|cn2rns2r|, which 
disappears at r = + n)K. The total number of gaps in the spectrum is three in the case 
r / (i + n)K, £ G (-00,^3,-] U [£2,-,£i,-] U [fi,+ ,f2,+] U [£^3,+ , 00), while for r = (i + n)K 
there are only two gaps, £ G (— oo,<?3^_] U [£2 -,£2,+] U [£^3^+, 00). According to (|3.15p . (|3.16p 
and (IS.lOp . the edges £i^e of the internal (i = 1, 2) and external {i = 3) allowed bands are 

fi,,(r) = ev^, f2,e(r) = e^/k'^ + e{r) , ^3,.(r) = e y^l + e(T) , (7.1) 
where e = it, and the eigenstates have a form ^i ,:{x,T) = ee*'^?'('^^'i/'i(^-)^ > 

S'l^'i,,(x,T) = £i^e'^i^e{x,T). 

In the context of the physics of conducting polymers, for example, the internal bands are 
referred to as the lower, [£2-,£i~], and upper, [£"1^+, iS2,+], polaron bands; the upper external 
band, [£3^^,00), is called the conduction band; the lower external band, {-co, £3^^], is referred 
to as the valence band [31] . In general case for eigenstates ()3.19p we have 

Si^l,S^Ax,T) = £,{a,T)^l^Suei^,^), £e{a,T) = ey'E{a)+e{T), (7.2) 

where E(a) for internal and external bands is given by Eqs. ()2.5p and (12. 6p . 

Since HbcLG = Si does not distinguish index it of the wave functions within the allowed 
bands, each corresponding energy level is doubly degenerate. Six edge states for t / (^ + n)K 
are singlets. In the case r = (^ +n)K, four edge states with energies £ = and ±1 are singlets. 
Zero energy states ^'i^e form a doublet in this case, as it happens for any other energy level inside 
any allowed band. 

The described degeneration in the spectrum of Si indicates that the BdG system might possess 
its own nonlinear supersymmetric structure. This is indeed the case. First of all, from Table [T] 
we see that there are three operators, TZai, T(T2 and TZTcr^, which commute with 5*1, and square 
of each equals one. Hence, each of them may be chosen as a Z2 grading operator for the BdG 
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Figure 3: Spectrum of Hbiig = Si possesses symmetries <Se(a,r) = £^{a,—T) = £e{oi,T + K), 
<?e(a, + r) = <5e(a, — r), and <5_(q, t) = — (?+(a, r). Horizontal line shows a spectrum for 
some value ofr, < r < K. The allowed (forbidden) bands on it are presented by thick green 
(thin red) intervals, whose points are distinguished by the parameter a, see Eq. (|7.2p . Curves 
indicate the edges of the allowed bands (j7.ip . The point £e{Vi + iK', ^K) = corresponds to a 
doubly degenerate energy level in the allowed band [— /c',A;'], that is formed by the two merging 
at T = iK internal allowed bands. 



system. There are three more, nontrivial local integrals of motion for HsdG- One is the second 
order operator T-L. This, however, is not interesting from the viewpoint of a supersymmetric 
structure since it is just a shifted square of HsdC = Si, V. = Sf — ^(t). Then we have a third 
order integral Li = Ci, that has been identified before as the Lax operator for the self-isospectral 
Lame system 71. Finally, the fourth order operator Qi = SiCi is also identified as a local integral 
of motion. Note that the both integrals Ci and Qi distinguish the states inside the allowed bands, 
which differ in index it . On distinguishing the states with iS = to be present in the spectrum 
if r = (i + n)K, see a discussion at the end of the previous section. Further nontrivial, but 
nonlocal integrals may be obtained if we multiply local integrals by the operators TZai, Ta2 and 
TZTcr^. Then, as in the case of the self-isospectral Lame system, different choices for the grading 
operator lead to distinct identifications of Z2 parities of the integrals. 

For the sake of definiteness, let us choose F = TZai, and assume first that r 7^ (^ + n)K. Other 
two possibilities for the choice of F may be considered in an analogous way. If, additionally, we 
restrict our analysis by the integrals that do not include in their structure a nonlocal in r operator 
T, we get two Z2-even (commuting with F) integrals in addition to HbcIG = Si, namely, TZai 
and TZaiSi. The four Z2-odd (anticommuting with F) integrals are Ci, Qi, C2 = iTZaiCi and 
Q2 = ilZaiQi- All these integrals are hermitian operators. It is interesting to note that a nonlocal 
integral TZai Si is related to one of the diagonal nonlocal operators from (16. lip . IZuiSi = a^S. 
A nonlocal diagonal operator Q2 also may be related to ()6.1ip . Q2 = QS^ + C(r)S'(5f — £{t)). 
Since, however, integrals TZaiSi and Ga are just the integrals TZai and multiplied by the BdG 
Hamiltonian Si, we can forget them as well as Ti. We obtain then nontrivial (anti) commutation 
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relations of the nonlinear BdG superalgebra, 

[IZai , Ca] = -2ieabCb , Cb} = 25„bP(5i , r) . (7.3) 

Here, in correspondence with Eqs. ()5.19p . (|5.2ip and ()6.5p . P{Si,t) is the six order spectral 
polynomial of the BdG system, 

r) = {Si - e{r)){S! - e(r) - k'^){S! - e(r) - 1) , (7.4) 

whose six roots correspond to the energy levels (17. ID . 

Superalgebra ()7.3p has a structure similar to that of a hidden, bosonized supersymmetry [U] 
of the unextended Lame system (12. ip . which was revealed in [38]. There, the role of the grading 
operator is played by a reflection operator TZ, the matrix integrals Ca are substituted by the Lax 
operator —iV{x), see Eq. (15.15j) . and by 'R.V{x). The six order polynomial P{Si,t) of the BdG 
Hamiltonian Si is changed there for a third order spectral polynomial P{H), see Eq. (j5.19p . 

We have seen that the structure of the BdG spectrum changes significantly at r = (^ +n)K. 
Essential changes happen also in the superalgebraic structure. Indeed, from (|6.8p it follows that 
[Si, Q2] = 2ia^eabC{T)e{T) , i. e. in a generic case Q2 does not commute with HsdG- In contrary, 
for r = (^ + n)K this is an additional nontrivial, second order integral of motion of the BdG 
system. This integral, like the third order integral Li, also distinguishes the states marked by 

the index it inside the allowed bands, Q2^± Si e ~ ^V\J Pqyi^i^)) ^± Si o where rj is the same 

as in (|2TT]1 and dOOjl . i.e. ry = -1 for < £^ < k'"^ and 77 = +1 for > 1, while PqyiE) is a 
polynomial that appeared earlier in (j5.33p . i. e. PQy{E) = {E — k''^){E — 1). In this case, Li is 
not independent integral for the BdG system anymore since here Li = —S1Q2 in correspondence 
with (I5.36p . Integral Q2 anticommutes with TZai and TZTcr^- Let us choose, again, F = TZai, 
and denote Qi = Q2 and Q2 = iTQi. Instead of (|7.3p . we get a nonlinear superalgebra of the 
order four, 

[nai, Qa] = -2ieabQb , {Qa, Qb} = 25abPQ{Sl) , (7.5) 

where Pq{Si) = {Sf - k'^){Sf - 1). 

It is interesting to see what happens with the Witten index in the described unusual super- 
symmetry of the BdG system with the first order Hamiltonian. One can construct the eigenstates 
of the grading operator F = TZai, 

r¥'\x;a,T) = -e^^'\x;a,T), ^^'\x;a,T) = ^"^^^^^(x, r) + e^'^''("'^)^^^5^^,(x, r) . (7.6) 

For any energy value inside any allowed band (including iS = in the case of r = (i +n)K), we 
have two states with opposite eigenvalues of F, and these contribute zero into the Witten index 
Aw = TrF, where trace is taken over all the eigenstates of the grading operator F. On the other 
hand, the edge states ^'i^e(x,r) are singlets. They are also the eigenstates of F. The eigenstates 
of opposite energy signs have opposite eigenvalues, +1 and —1, of the grading operator. As a 
result, we conclude that the Witten index Aw in such a supersymmetric system equals zero for 
any value of r [i.e., for t 7^ (| + n)K when there are no zero energy states in the spectrum, and 
for T = (| + n)K when the spectrum contains a doublet of zero energy states], like this happens in 
the self-isospectral Lame system with the second order supersymmetric Hamiltonian. The same 
result Aw = is obtained for the choices F = Tcr2 and F = IZTa^. 

Finally, it is worth to notice that in accordance with the structure of superalgebra (|7.3p . the 
third order matrix BdG supercharges Ca annihilate all the six edge eigenstates of Hsdc = 5*1 
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in the case of r 7^ + n)K. In special cases r = (2 + n)K. a central gap disappears in the 
spectrum, and, consistently with (|7.5|) . all the remaining four edge states are the zero modes of 
the second order matrix BdG supercharges Qa- In other words, the spectral changes that happen 
in the BdG system at special values of the parameter t = + n)K, which correspond to a zero 
value of the bare mass thq in the GN model (jl.2p . are reflected coherently by the changes in its 
super algebraic structure. 



8 Infinite period limit 

Let us discuss now the infinite period limit of our self-isospectral Lame and the associated BdG 
systems, i. e. the case when the period 2K tends to infinity. 

K — )• 00 assumes El /c — t- 1, /c' — ?■ 0, K' — ?■ ^tt, and relations (jA.Sp . and (jB.SP have to be 
employed. According to (jB.Sp and (|B.9p . a limit for a quotient of functions is also well defined, 

@{u) cosh(u) 

^"1?H = u,v€C. (8.1) 

fc^-i ld>[v) cosh(uj 

Periodic Lame Hamiltonian (|2.ip transforms in this limit into a reflectionless one-gap Poschl- 
Teller Hamiltonian 



When the limit K — )• c« is applied to the self-isospectral system (j3.8p . we assume that a shift 
parameter r remains to be finite. As a result we get a self-isospectral non-periodic FT system, 

TipTix) = diag {Hr{x),H^r{x)) , (8.3) 

where Ht{x) = Hpt{x + r) and H-t-{x) = HpTix — r). In what follows we trace out how the 
peculiar supersymmetry of the self-isospectral Lame system transforms in the infinite period limit 
into the super symmetric structure of the system (18. Sp . which was studied recently in [40j. 

Since the super-partners in (j8.3p are the two mutually shifted copies of the same PT system, 
it is clear that the limit does not change the Witten index : it remains to be equal zero as in the 
periodic case. In general, however, the index may or may not change depending on the concrete 
form of the self-isospectral Lame system to which the limit is applied. For instance, in the case 
of the system with superpartners H{x) and H{x + K) [see a remark just below Eq. ()3.7p ]. the 
infinite-period limit gives, instead of ()8.3p . a supersymmetric system with one superpartner to 
be the FT system ()8.2p . while another one (which is a limit of H{x + K)) to be a free particle 
Hq = —-^ + 1. Superpartner potentials in such a supersymmetric (but not self-isospectral) 
system are distinct. The only difference of the spectrum for the system ()8.2p from that for Hq 
consists in the presence of a unique bound state, see below. Consequently, the Witten index 
changes in the infinite period limit, by taking a value of the modulus one. If in the system p.8p 
one takes r = t(K) such that r — t- 00 for K — t- 00, the limit produces then a trivial self-isospectral 
system composed from the two copies of the free particle Hamiltonian Hq. In such a case, the 
Witten index does not change in agreement with (18. 3p and (18. 2p . 

The listed examples also mean that the shifts for the period, in a sense, 'interfere' with the 
infinite period limit. Self-isospectral Lame system composed from H{x+) and H[x-) is equivalent, 
for instance, to a system with super-partner Hamiltonians i?(x+) and H{x- -\- 2K) 0. If before 



^■^Any of these four limits assumes three others. 



The second system, however, is characterized by another phase (|3.14|) with r changed for r — K. 
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taking a limit we do not 'eliminate' the period 2K shift in the second subsystem, we will obtain 
a (not self-isospectral) system with super-partners H-^ and Hq instead of ()8.3p . 

Let us return to the symmetric case of the self-isospectral Lame system (jS.Sp , whose infinite 
period limit corresponds to the self-isospectral PT system (jS.Sp . All the energy values (j2.5p of 
the valence band transform into zero in the infinite period limit because of k' — )• 0, i. e. all 
this band shrinks just into a one energy level E = Q for the system (18. 2p . In conformity with 
this, all the Bloch states ()2.3p of this band, including the edge states dnx and cnx, turn into 
a unique bound state ^^^-^ ^ of S = for PT system El. Then the states 1/ cosh(x it r) form 
a supersymmetric doublet of the ground states for self-isospectral system (j8.3p . The doublet 
of the edge states sn {x ± r) of the system (13. 8j) transforms into a doublet of the lowest states 
tanh(x±r) of the energy £^ = 1 in the scattering sector of the spectrum for (|8.3p . It is interesting 
to see how the eigenstates with E > \ m the scattering sector of the PT system originate from 
the Bloch states (|2.3p . The energy (|2.6p as a function of the parameter /?, which in the infinite 
period limit takes values in the interval < /? < ^, reduces to E{ij3) = ^^^^ ^ > 1. The states 

()2.3p transform into = cos/3 (tanhx it ztan/3)exp(=pzxtan/3). Denoting tan/3 = k > 0, 

we obtain E = 1 -|- k^, and the states ^'^(x) take the form of the scattering eigenstates of the 
PT system, — > ^^^{x) = — - tanhx)e=^*^. 



We have 



cosh x+ 

for function ()3.2p . cf. Eq. (5.17) in [40] . In correspondence with (j3.4p . this is a nonphysical 
eigenstate of Hj^ of eigenvalue —1/ sinh^ 2r. Function A(x; r) in the form (I4.ip transforms into 

A(x;r) > Ar{x) = coth 2t + tanh x_ — tanhx^, (8.5) 

while Eq. ()4.2p gives, equivalently, 

k^l 2 

A(x;t) > At-(x) = — — — h tanh 2r tanh tanh X-L . (8-6) 

sinh 4r 

Non-periodic superpotential (gap function) (18. 5p corresponds to the DHN kink-antikink baryons 
[2]. For the first order intertwining operator we have 

V{x;T)^^-AAx)^Xr, (8.7) 
cf. (2.26) in [40]. It is the operator that appears in the limit structure of the supercharges Sa, 

Si ( ^ I = SpTi , 5*2 SpT.2 = icTsSpT.i ■ (8.8) 



Xr 

For the second order intertwining operator (j5.6p . 

A{x; t) ^ A^rAl = Yr , (8.9) 

where limK-s.oo T^{x + t + ^K; ^K) = limK^.oo T^{x + ^K; -r + |K) = ^ - tanh x+ = At-{x), 
and A-T is obtained via the change r — )• — r. A limit of the second order integrals ()5.7p is 



Ql > [ y J" ) = Q^'T,! , (52 > QPT,2 = i(T3QpT,l , (8.10) 



^^The states (|2.3p for the valence band should be 'renormalized' (divided) by a constant 0(K)/0(O) to cancel 
the multiplicative factor that diverges in the limit K — >■ oo in correspondence with (|8.1[) . 
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cf. Eq. (2.18) in [30]. The first order operators Aj- and A^r factorize also the self-isospectral pair 
of the PT Hamiltonians, = A^- At, H^r = A^rA^-r^ 

as well as a free particle Hamiltonian, 

= aXa^ = A^_^A_r- 

The phases that appear in the action of the intertwining operators r) and A{x\ r) on the 
super-partner's eigenstates, see Eqs. (|3.1ip and (jS.lOp . transform into 

^iHP{c.,r) g-2*r . "^k " COth 2r ^ gi^^(a,r) g-2ikr _ ^^ -y^^ 

s2 + coth^ 2r ' 



They are associated with the action of the intertwining operators and 1^ on the eigenstates 
of super-partner systems Hr and -ff-r, and appear in the structure of the eigenstates of the first, 
(jS.Sp . and the second, (jS.lOp . order integrals of the self-isospectral PT system [40] . 

By employing a relation 2^(2:4.) = T)'^ {x\t^A(x\t^ — r)P(x; r) that follows from Eq. 

([STi]) . we find that 



V^x^)"^ A^^Al^Z^, (8.12) 
cf. (2.24) in [40]. For the limit of the Lax integrals we get then 

Li ~^ Q Z ) ^ LPT,1 , L2 -^^^ LpT,2 = o-sLpT,! ■ (8.13) 

Finally, for a constant C(r) = ns2Tnc2rdn2r that appears in the superalgebraic 
(anti) commutation relations of our system we obtain 

C(t) COth 2r = C2r , (8.14) 

cf. the first term in Eq. (18. 5p . 

With the described infinite period limit relations, we find a correspondence between the 
supersymmetric structures in the self-isospectral one-gap Lame and PT systems. Particularly, 
applying the infinite period limit to the superalgebraic relations of the self-isospectral Lame 
system and making use of the described correspondence, one may reproduce immediately the 
superalgebraic relations for the self-isospectral PT system (jS.Sp . 

The same r-dependent constant C2T = coth 2r shows up in representation for superpotential 
(jS.Sp and in the superalgebraic structure for the self-isospectral non-periodic PT system ()8.3p 
due to relation (j8.14p . Notice, however, that corresponding functions of a shift parameter, z(t) 
and C(t), which appear in the periodic system, are different. In the next section we will return 
to this observation. 

The infinite-period limit of the second order intertwining operator 3^(2;; r) may be found by 
employing relation ()5.6p . 

lim y{x; r) = -Y, - Car^r • (8.15) 

K—>-oo 



It plays no special role in the supersymmetric structure of the self-isospectral PT system (18. Sp . Let 
us, however, shift x — )• x — r in (jS.lSp and then take a limit r — )• 00. Such a double limit procedure 
applied to the self-isospectral Lame system T-L produces a non-periodic supersymmetric system 
v. = diag {Hpt{x),Hq{x)) that is composed from the PT system ()8.2p and the free particle 
Hq = —-^ + 1. Operator 3^(x;r) in such a limit transforms into the second order operator 
y(^) ~ lx('^ ^^^^ intertwines HpT with Hq, y{x)HpT{x) = HQ{x)y{x). The kernel of y 

is formed by singlet eigenstates 1/ coshx {E = 0) and tanhx {E = 1) of the PT system Hpt{x)^ 
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cf. the discussion of the kernel of y{x; ^K) in section [5j Hermitian conjugate operator y^{x) 
intertwines as {x)Hq{x) = HpT{x)y^ (x), and annihilates the eigenstate 1 of the lowest energy 
E = 1 and a non-physical state sinhx of zero energy in the spectrum of Hq. Integrals Sa, Qa 
and La transform in such a double limit into the integrals of the super symmetric system 'H, 



and 5*2 — 5- S2 = icrs-h, Q2 ~^ Q2 ~ ^'^39i> L2 ^ h = crsh, where Aq = limT-_^oo ^t(2; — t) = 



Non-periodic superpotential (gap function) A(x) = tanhx that appears in the structure of the 
first and second order intertwining operators as well as in that of the integrals (j8.16p corresponds 
to the famous CCGZ kink solution [ig [30] of the GN model. 

Prom the total number of seven integrals of motion (j6.ip and a^, each of which can be used 
as a grading operator for self-isospectral Lame and PT systems, in the described double limit 
survive only three: in addition to the obvious operator 173, nonlocal operators TZ and TZcr^, are 
also the integrals for super symmetric system Ti. The last two operators originate in the double 
limit from the integrals TZT and TZTc^. Having in mind this correspondence, Table [1] still may 
be used for identification of Z2 parities of the integrals Sa, qH and la, and it is not difficult to 
obtain corresponding forms for superalgebra for each of the three possible choices of the grading 
operator in this case, see [39t HH], 

Let us look what happens here with the Witten index. As we discussed at the beginning 
of this section, the only asymmetry between the spectra of the superpartner Hamiltonians HpT 
and Hq is the presence of the zero energy bound state in the first super-partner system, which 
is described by the eigenstate (1/ coshx,0)-^ of the super symmetric system T-i. The doublet with 
E = 1 is formed by the eigenstates (tanhx,0)^ and (0,1)^. The first state is an eigenstate of 
all the three operators C73, TZ and TZa^ with the same eigenvalue +1, while for the second and 
third states the eigenvalues are, respectively, +1, —1, —1 and —1, +1, —1. All the forth-fold 
degenerate energy levels in the scattering part of the spectrum with E > 1 contribute zero into 
the Witten index Avi^ = TrP. As a result, for all the three choices of the grading operator for 
non-periodic supersymmetric system ?i we have consistently |Avk| = 10- 

On the other hand, the first order matrix operator si is identified here as a limit of the 
BdG Hamiltonian Hsdc = Si. As may be checked directly, operator IZa^ commutes with si in 
accordance with Table [1] if to take into account the correspondence between nonlocal integrals 
discussed above. Therefore, it can be identified as a grading operator for a peculiar supersym- 
metry of the BdG system with the Hamiltonian h^dG = ■^i, in which the second order integral 

and the nonlocal operator ilZa^q^ are identified as the odd supercharges, and li = —siq^, 
cf. (j5.36p . Corresponding superalgebra has a form (|7.5p with obvious substitutions. The state 
(1/ coshx,0)^, is a unique zero mode of the first order matrix hamiltonian si, while two states 
(tanha;, ±1)-^ are the singlet eigenstates of si of the eigenvalues ±1, which are also the eigenstates 
of the grading operator IZa^ of the eigenvalue —1. 

Thus, the modulus of the Witten index changes from zero to one for the supersymmetries of 
the both, second, and first, hBdG = ^i; order systems. This refiects effectively the changes 

^^Avi/ takes values +1 for F = aa and TZ, and —1 for TZa-i. A difference in sign is not important, fiowever, since 
it can be removed by changing a sign in definition of the grading operator in the last case. 




(8.16) 




0, 
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in the spectrum that happen in the described infinite-period hmit of the self-isospectral second 
order Lame and the associated first order BdG systems. 

9 Extended supersymmetric picture and Darboux dressing 

Let us discuss now another interesting aspect of our self-isospectral periodic supersymmetric 
system in the light of the infinite period limit. As it was shown in [40], the supersymmetric 
structure of the non-periodic self-isospectral system (|8.3p has a peculiar property : all its integrals 
can be treated as a Darboux-dressed form of the integrals of a free particle system Hq{x). We 
clarify now what corresponds here, in the periodic case, to the Darboux-dressing structure of 
the self-isospectral PT system (|8.3p . For that, we extend a picture related to the intertwining 
operators and the Darboux displacements associated with them. 

Consider along with our self-isospectral supersymmetric Lame system (j3.8p . ^{{x) = 
diag {H{x + t),H{x — r)), its copy shifted for the half period, ^{{x + K) = diag [H{x + K + 
t),H(x + K — r). Any two of the four (single-component) Hamiltonians may be connected by 
intertwining relation of the form ^j)H{^ + /i) = H{£^ — /i). Putting ^ = x + ^(n + r2) 

and = ^(ri — t"! 7^ ^2 + 2Kn, we present this relation in a more appropriate form 

V{x + i(ri + r2); \{ti - T2))H{x + n) = H{x + T2)V{x + ^{n + ra); i(ri - ra)) . (9.1) 

Here ri and T2 take values in the set {— r, r, — r + K, r + K}, and supersymmetric Hamiltonians 
n{x) and n{x + K) may be related by VTiix + K) = 'H{x)V, V'^n{x) = n{x + K)pt, where 

V = diag {V{x + r + iK; ^K), P(x - r + ^K; ^K)) . (9.2) 

In general case, if any two Hamiltonians h and h are related by intertwining operators D and 
L>t, Dh = hD, hD'f = D^h, and if J is an integral for h, [h, J] = 0, then the operator J = DJD^ 
is an integral for h. The system ^(x) is characterized by the set of local integrals of motion 
J{x) = {(T3, Sa{x), Qa{x), La(x)}, while the system 'H{x + K), is described by the same but 
shifted set, J{x + K). Identifying '}i{x + K), 'H{x) and V with h, h and D, respectively, we 
find that J = T)J[x + K)P^ = J{x)'H{x). In other words, the Darboux dressed integral of one 
system is just the corresponding integral of another, displaced self-isospectral periodic system, 
multiplied by its Hamiltonian. Nonlocal operators (16. ip . which are the integrals for 1-L{x), are also 
the integrals of motion for the displaced system 'H(a; + K). Then one finds that a similar relation 
is valid also for these nonlocal integrals as well as for nontrivial diagonal nonlocal integrals ()6.1ip . 
The only difference is that for all the integrals that contain a factor TZ, including (16. lip , there 
appears a minus sign, like in 'DS{x + K.)!)^ = —S{x)'H{x). Notice also that the Darboux dressed 
form of the trivial integral 1 (that is a unit two-by-two matrix) for the displaced system T-L{x + K.) 
coincides with the Hamiltonian Hix), VIV^ = 'H{x). 

Since the both self-isospectral supersymmetric systems are just two copies of the same periodic 
system shifted mutually in the half period, the described picture is not so unexpected. Let us 
look, however, at this result from another viewpoint. In the infinite period limit, supersymmetric 
systems 'H{x) and 'H(x + K) transform, respectively, into (jS.Sp and 

1^0 = diag (i/0,^0), (9.3) 

where Hq = —-^ -|- 1 is a (shifted for a constant additive term) free particle Hamiltonian. In 
other words, the infinite period limit of the system T-L{x + K) is given by the two copies of the 
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free non-relativistic particle. As we have seen, the infinite period hmit apphed to the integrals 
of the self-isospectral system T-L^x) produces corresponding integrals of the self-isospectral PT 
system (j8.3p . The infinite period limit of the integrals of the system ^{{x + K) may easily be 
obtained just by taking a limit x — )• oo of the integrals of the self-isospectral PT system (18. 3p . 
For nontrivial local integrals we find 

+ K) — > -i-^(T2 - C2rO-l = si , S'2 (x + K) — ^ §2 = ^CTsSi , (9.4) 
ax 

Qa{x + K) ^ (-1)°+V„, • ^0 , ^1(2; + K) ^ -i^- -^0 = 11, L2{x + K) ^ ^2 = cy^h ■ (9-5) 

ax 

The obtained operators are the integrals of motion for the trivial free particle super symmetric 
system (j9.3|) . They correspond to the obvious integrals a a-, and to the products of them with —i-^ 
and 1-Lq. System (j9.3p is intertwined with the self-isospectral PT system (jS.Sp by the infinite period 
limit of the operator V diag (A^,^_^) = Poo, i^oo^o = 'HptD^, ?^o^L = DtoHpT. 

If Jo is some integral for T-Lq, then DoqJqHqD]^ = DooJqDIo'Hpt- Taking into account (j9.4p 
and ()9.5p . the nontrivial local integrals SpT,a, QpT,a and LpT,a of the self-isospectral PT system 
(jS.Sp may be treated as a Darboux dressed form of the integrals for the free particle system TIq, 
namely, of Sa, (Ta and —iZa-^, where Xi = 1 and X2 = cs- 

It is interesting to note that the first order integral of for instance, si, may also be treated 
as a Hamiltonian of a free relativistic Dirac particle of mass C2t- Then its Darboux dressed form 
is a non-periodic BdG Hamiltonian 

SpT,i = -i-i-cr2 - Ar{x)ai, (9.6) 
ax 

see Eqs. (jS.Sp and ()8.5p . Comparing (j9.6p with the structure of si in ()9.4p . we see that a gap 
function At-{x) is effectively a Darboux dressed form of a free Dirac particle's mass C2t- The 
periodic BdG Hamiltonian HsdC = "5*1 may be treated then as a periodized form of (j9.6p . like the 
Lame Hamiltonian may be considered as a periodized form of the PT Hamiltonian, see [31]. It is 
worth to stress, however, that a reconstruction of a crystal structure on the basis of a non-periodic 
kink-antikink system is not direct and free of ambiguities : in the previous section we already 
noted that two different basic functions of the shift parameter in the self-isospectral Lame and 
associated BdG systems correspond to the same function in the non-periodic case. 

Another interesting observation can be made on a genesis of the non-local integrals ()6.1ip . 
For self-isospectral Lame and PT systems, the reflection operator TZ and da, a = 1,2, are 
not integrals of motion, but the product of any two of these three operators is an integral 
of motion. For supersymmetric free particle system (j9.3p . however, each of these three op- 
erators is an integral of motion. One finds then that the infinite period limit of the in- 
tegral (TsQ, — ^ diag {JZYt^IZY-t^) = a^QpT is exactly a Darboux dressed form of the 
reflection operator TZ, DooTZDIo = ct^Qpt- Or, alternatively, an integral Qpt for the self- 
isospectral PT system is a dressed form of the nonlocal diagonal integral TZa^. An analogous 
relation exists also for the inflnite period limit of another nonlocal diagonal integral from (16. lip . 
Doo{—iTZa2Si)Dlo = Spt ■ T-Lpx, where Spx = diag {TZXr,TZX-r)- 

We conclude that the described Darboux dressing structure of the self-isospectral PT system, 
observed earlier in HQ], originates from, and is explained by the properties of the self-isospectral 
periodic one-gap Lame system. 
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10 Discussion and outlook 



To conclude, let us discuss the obtained results from the physics perspective and potential appli- 
cations and generalizations. 

Usual supersymmetric structure of the kink-antikink as well as of the kink crystalline phases 
of the GN model is known for about twenty years. However, such a structure with the first order 
supercharges and Z2 grading provided by the diagonal Pauli matrix does not explain or reflect a 
peculiar, finite-gap nature of the corresponding solutions. It does not reflect either a restoration 
of the discrete chiral symmetry at zero value of the bare mass in the GN model, when the kink- 
antikink crystalline condensate transforms into the kink crystal. The both aspects are explained 
by the exotic nonlinear supersymmetric structure we revealed here. The finite-gap nature is 
reflected by the Lax integral incorporated into a nonlinear supersymmetric structure alongside 
with the first and second order supercharges. A restoration of the discrete chiral symmetry, on 
the other hand, is reflected by structural changes that happen in nonlinear supersymmetry at 
the half period shift of Lame superpartner systems, when a central gap in the spectrum of the 
associated BdG system disappears. We showed that the first order BdG system has its own 
supersymmetry, which can be revealed only with the help of the nonlocal grading operators inves- 
tigated in Section 6. The disappearance of the middle gap in the BdG spectrum is accompanied 
by emergence of the new, nontrivial second order integral of motion in the first order system 
(while the BdG Hamiltonian has no such integral in the kink-antikink crystalline phase). 

The aspects related to the infinite period limit we investigated in sections 8 and 9 may be 
useful for understanding of some puzzles related to a computation of the Witten index in some 
supersymmetric field theories when a system is put in a periodized box [49]. 

Recently, perfect Klein tunneling in carbon nanostructures was explained in [50] by unusual 
supersymmetric structure with the first order matrix Hamiltonian. We believe that the super- 
symmetry we investigated here, particularly in Section 7, may also be useful in the study of other 
phenomena in graphene, where the dynamics of charges is governed by the effective first order 
Dirac Hamiltonian. 

It would be interesting to clarify whether the twisted kink crystal of the GN model with 
continuous chiral symmetry, that was found in [7', 'H], could be obtained by supersymmetric 
constructions similar to those from section 3. 

We treated A that appears in the structure of the second order intertwining operator B{x; r. A) 
of a general form (j5.3p as a kind of a virtual shift parameter. One could extend the picture by 
reinterpreting Eqs. ()5.ip and ()5.2p as intertwining relations for three Lame systems, H(x + ti), 
H{x + T2) and H{x + T3), where ti = t, T2 = t + 2X and T3 = — r. Then we would get an 
extended self-isospectral system of three super-partner Lame Hamiltonians. Employing relation 
of a form ()9.ip . one could further extend the picture to obtain a self-isospectral system with 
n > 3 superpartners H{x + ti), . . . , H{x + Tn)- When the shift parameters are such that r„ = ti, 
the corresponding intertwining operator of order n would reduce to an integral for the system 
H{x + Ti). It is in such a way we identified, in fact, the third order Lax operator V{x + r) 
for the system H(x + r). The interesting questions that arise are then: what is a complete set 
of integrals and what kind of supersymmetry we get for such an n-component self-isospectral 
system? Particularly, what is the nature of the above-mentioned integral of motion of the order 
n for n > 3? What is a relation of such extended supersymmetric systems with the GN model 
and what physics could be associated with them? 

^^It is this first order system that really describes the corresponding crystalline phases in the GN model while 
the second order Lame system is related to it as the Klein-Gordon equation is related to the Dirac equation. 
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Appendix A: Jacobi elliptic functions 

We summarize here some properties and relations for Jacobi elliptic and related functions. For 
details, see, e. g., [29l [5T] . 

In notations for these functions we suppress a dependence on a modular parameter < A: < 1, 
snx = sn{x\k), etc., when this does not lead to ambiguities. On the other hand, a dependence on 
a complementary modulus parameter < k' < 1, k' = {1 — A:^)^/^, is indicated explicitly. We use 
Glaisher's notation for inverse quantities and quotients of Jacobi elliptic functions, ndx = 1/dnx, 
nsx = 1/snx, ncx = 1/cnx, sex = snx/cnx, etc. 

The basic Jacobi elliptic functions are the doubly-periodic meromorphic functions snu, cnu 
and dnu, whose periods are (4K, 2iK'), (4K, 2K + 2iK') and (2K, 4iK'), respectively, snu is an 
odd function, while cnn and dnu are even functions, which are related by identities sn^n+cn^n = 
1, dn^u + k'^sn^u = 1, k'^cn?u + k'"^ = dn^n, /c'^sn^n + cv?u = dn^ti, and whose derivatives are 
-^snu = cntidnti, ^cnn = —snu dnu, -^dnu = —k'^snucnu. They have simple zeros and 
poles at 

snn:0,2K; cnu:K,-K; dnn : K + iK', K - iK' , (A.l) 
snu, cnu : iK', 2K + 2iK'; dnu : iK' , -iK' , (A.2) 
respectively, modulo periods. Here 

K = K(A;) = / (A.3) 

^ ' Jo ^(l-x2)(l-A;2x2) ^ ^ 

is a complete elliptic integral of the first kind, and K' = K(fe') is a complementary integral, which 
are monotonic functions of k in the interval < A: < 1 : dK/dk > 0, dK! /dk < 0. In the limit 
cases A; = and k = 1, elliptic functions transform into simply-periodic functions in a complex 
plane, 

A; = 0, A;' = 1 : K = ^tt, K' = oo , snu = sinu, cnu = cosu, dnu = 1, (A. 4) 
A; = 1, A:' = : K = oo, K' = ^vr , sntt = tanhu, cnu = dnu = \ . (A. 5) 



The addition formulae are 

S-j- — (suCydu -\- s^c^d^), C-|_ — i^u^v s^iS-yd-^d^) , d-|^ — (d-^^d-^; k SuSyCuCu\ (A. 6) 

/X /X /i 

where s+ = sn (u + v), Su = snu, s^ = snu, c_(_ = cn (u + u), d+ = dn (u + v), etc., and 
/i = 1 — k^sn^u sn^v . The Jacobi's imaginary transformation is 

sn(iu|A;) = isn(u| A;')nc(u|A;'), cn(iu|A;) = nc(u|A;'), dn(iu|A;) = dn(u| A;')nc(u| A;') . (A. 7) 

From addition formulae and ()A.7p . one finds some displacement properties of Jacobi elliptic 
functions shown in Table [2j 
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Table 2: Displacement properties of Jacobi elliptic functions 



n 


u + K 


u + iK' 


u + K + iK' 


U + 2K 


u + 2iK' 


u + 2(K + iK') 


snu 


en undu 


^ nsu 


^ dnuncu 


—snu 


snu 


—snu 


cnu 


—k'snundu 


—i-^ dnunsu 


■ k' 

— new 


—cnu 


—cnu 


cnu 


dnu 


k'ndu 


—icn unsu 


ik'sn uncu 


dnu 


—dnu 


—dnu 



Appendix B: Jacobi Zeta, Theta and Eta functions 

The complete elliptic integral of the second kind is defined by 







E = E(fc)= / \ - ^dx. (B.l) 



1 — X 



It is a monotonically decreasing function, dE/dk < 0. The complete elliptic integrals K = K(/c) 
and E = E(A;) satisfy the first order differential equations ^ = ^^fr^, ^ = bom which 

an inequality k''^ < E/K < 1 and the Legendre's relation EK' + E'K — KK' = ^tt may be 
deduced, where E' = E(A;') is a complementary integral of the second kind. 
The incomplete elliptic integral of the second kind is defined as 

E{u) = / dn^udu, (B.2) 
Jo 

in terms of which E = E(K). This is an odd analytic function of n, regular save for simple poles of 
residue +1 at the points 2nK+{2m+l)iK' . Function E(n) is not an elliptic function. It possesses 
the properties of pseudo-periodicity, E(w + 2K) - E(u) = E(2K) = 2E, E(m + 2iK') - E(m) = 
E(2zK'), where in the first relation the second equality is obtained by putting u = — K. 
In terms of E('u), a simply-periodic Jacobi Zeta function is defined, 

Z{u) = E{u)-^u, (B.3) 

which satisfies relations ^^j^ = dn^-u — ^ , and 

Z(u + 2K) = Z(«), Z(u + 2iK') = Z{u) - if, Z{-u) = -Z{u), Z(K - u) = -Z(K + u),{BA) 

Z(0) = Z(K) = 0, Z(K + iK') = -i^. (B.5) 

Zeta function satisfies an addition formula 

Z{u + v) = Z{u) + Z{v) - k\nusnvsn{u + v) , (B.6) 

and obeys Jacobi's imaginary transformation 

ITU 

iZ{iu\k) = Z{u\k') + — — - dn{u\k')sc{u\k') , (B.7) 
2xCIC 

from which one finds Z[u + iK!) = Z(?i)-|-nsncn'udn'u For the limit values of the modular 

parameter, /c = and = 1, we have 

Z(u|0) = 0, Z('u|l) =tanh'u. (B.8) 
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In terms of Z(n) = Z{u\k), the Jacobi Theta function @{u\k) is defined as 



e(n) = 6(0) exp 



Z{u) du 



(B.9) 



This is an even, Q{—u) = 0(n), integral periodic function of period 2K, whose only zeros are 
simple ones at the points of the set 2nK + (2m + l)iK'. It satisfies a relation Q{u + 2iK') = 
— i exp {—i^u) where q = exp(— vrK'/K). Notice that sometimes Jacobi's Theta function 

is defined by the Fourier series, 



e{u\k) = i)^{v) , Mz) = 1 + 2 ^(-l)"g"' cos(2 



nz) 



n=l 



-KU 

2K 



Then Z function can be defined by the logarithmic derivative, 

Z(n) = Alne(n). 
du 



(B.IO) 



(B.ll) 



In correspondence with definition (jB.lOp . a constant in ()B.9p is fixed as 0(0) 
Jacobi Eta function H(m) is defined in terms of the Theta function. 



B{u) = -iq^/^ exp (i^) e(n + iK') 



2Kfc' 



(B.12) 



This is an odd, II(— u) = — H(m), integral periodic function of period 4K, which possesses simple 
zeros at the points of the set 2nK + 2miK'. Some properties of the Eta and Theta functions are 
summarized in TableO where M{u) = exp N{u) = exp {—i^) q~^ ■ For particular 

Table 3: Parity and some displacement properties of Jacobi and H functions 



u 


—u 


U + 2K 


u + iK' 


u + 2iK' 


u + K + iK' 


U + 2K + 2iK' 


e(n) 


e(n) 


eiu) 


iM{u)R{u) 


-N{u)e{u) 


M(n)H(n + K) 


-N{u)e{u) 


H(n) 


-H(n) 


-R{u) 


iM{u)eiu) 


-N{u)R(u) 


M(n)e(u + K) 


N{u)R{u) 



values of the argument we also have H'(0) = 2KH(K)e(0)e(K), e(K) = H(K) = 

Jacobi Theta function satisfies a kind of addition theorem, 

e{u + v)Q{u - v)e\0) = Q\u)Q\v) - h2(u)h2(u) . (B.13) 

The basic Jacobi elliptic functions may be represented in terms of and H functions. 



H(n) 6(0) H(n + K) 6(0) 

snu = „ , , • — — ^ , cn n — 



6(n) H'(0) ' 



e{u) H(K) 



@{u + K) 6(0) 
^^^ = ^6(^-6M- ^""-''^ 



Under the complex conjugation, all the Jacobi elliptic functions as well as H, 6 and Z satisfy 
a relation (/(z))* = /(z*). 
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